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■ ABSTRACT 

o : 

. Context. 

' Aims. We determined the K s band luminosity function (LF), and inferred the corresponding stellar mass function, of cluster galaxies at redshift 
, z — 1.2, using near-infrared images of three X-ray luminous clusters at z = 1.11, 1.24, 1.27. 
O |' Methods. The composite LF was derived down to M*+4, by means of statistical background subtraction, and is well described by a Schechter 
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function with A** = 20.5+j 4 and a = -LO+jJ-j. Using available X-ray mass profiles we determined the MIL ratios of these three clusters, which 
tend to be lower than those measured in the local universe. Finally, from the K s band composite LF we derived the stellar mass function of 
cluster galaxies. 

j-j Results. With these data, no significant difference can be seen between the cluster galaxies LF and the LF of field galaxies at similar redshift. 
• • ■ We also found no significant evolution out to z — 1.2 in the bright (< M*+4) part of the LF probed in this study, apart from a brightening of 
^ = 1.3 mag of the characteristic magnitude of the high redshift LF. We confirm, and extend to higher redshift, the result from previous work that 
K> the redshift evolution of the characteristic magnitude M* is consistent with passive evolution of a stellar population formed at z > 2. 
t I ■ Conclusions. The results obtained in this work support and extend previous findings that most of the stars in bright galaxies were formed at 
' high redshift, and that A",-bright (M > 10" M Q ) galaxies were already in place at z — 1.2, at least in the central regions of X-ray luminous 
clusters. Together with recent results on the field galaxies stellar mass function, this implies that most of the stellar mass is already assembled 
in massive galaxies by z ^ 1, both in low and high density environments. 

Key words, galaxies: clusters: individual: RDCS J1252-2927, RDCS J0910+5422, RX J0848+4453 - galaxies: evolution - galaxies: formation 
-galaxies: luminosity function, mass function - cosmology: observations 
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1. Introduction 

Galaxy clusters are rare systems forming in the highest density 
peaks of large scale structure. In these special regions galaxy 
formation and evolutionary processes are expected to be faster 
with respect to the low density fields, thus making galaxy clus- 
ters a biased environment. On the other hand, clusters of galax- 
ies, particularly at high redshift (z - 1), provide a very conve- 
nient place for studying the evolution of massive galaxies. Not 
only do they contain high numbers of such objects, but these 
objects turn out to be so evolved (already at z - 1) that they 
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show a colour-magnitude sequence as clear as at lower red- 
shifts. Thus evolved galaxies in distant clusters can be easily 
identified even without complete spectroscopic follow-up. 

The study of massive galaxies has a relevant role in con- 
straining galaxy formation and evolution models, as different 
models provide different predictions for their assembly (in par- 
ticular in the redshift range [0 -f 1]). They could have rapidly 
formed their stars at high redshift and at the same time assem- 
bled their stellar mass, and then simply evolved passively as 
their stars aged. Alternatively, massive galaxies could have as- 
sembled on a longer time scale in a process of continuous merg- 
ing of smaller units until redshift < 1 . Comparison of these dif- 
ferent scenarios has proven to be a difficult task: even if merg- 
ing galaxies are observed, the relevance of the merging process 
in galaxy evolution and especially the epoch at which major 
mergers occur is still debated. 

Colours and spectra of massive galaxies at z — 1 show that 

there is a significant population of such systems already hosting 

mainly evolved stellar populations, both in the field (see below) 
j 1 1 j i I ' ' \ 

and in clusters (Stanford et al. ( 1997 1998), van Dokkum et al. 
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Jl999h. 


De Propris et al. 


il999), 


Rosatietal. (1999), . 


Stanford et al. (2002), Blakeslee et al. 


(2003), van Dokkum & Stanford (2003), 


Kodama et al. (2004), 


Lidman et alJ 




2004). iDe Lucia et alJ 


2004), 


Holden et alJ 



of evolved stellar populations is generally inferred from 
fundamental plane or colour-magnitude sequence evolution. 
These studies indicate that most of the stars in massive galaxies 
were formed at z > 2. At the same time several of them point 
out that the less massive the galaxy is, the more likely is the 
presence of a younger component in its stel lar population - the 
so called "downsizing" JCowie et alJ [T996) in galaxies hosting 
star formation. 

It is not possible, however, on the basis of spectrophotomet- 
ric analysis only, to rule out the possibility that these galaxies 
formed via merging of smaller galaxies with already evolved 
stellar populations even a short time before being observed. 
For instance, a passively evolving zero-point of the colour- 
magnitude relation does not imply that the galaxies formed 
long ago, but that the stars in the galaxies formed at high red- 
shift, possibly in smaller progenitors. 

In other words, while the underlying stellar populations 
can place constraints on the details of the star-formation his- 
tory, they cannot tell when a galaxy assembled. In fact, even if 
these massive galaxies appear to be passively evolving, several 

studies have noted that to some extent they can still be form- 

J I 1 

ing (or recently have formed) stars (for instance Nakata et al. 
ll200ll).lvan Dokkum & Stanford! (Eool.lHolden et alJ ll2005ct) . 

al 

and lTran et alJ ll2005l) ) or anyway subsequent episodes of star 
formation. Such secondary episodes of star formation are prob- 
ably correlated to cluster-related processes (accretion of field 
galaxies or groups and cluster merging), since galaxies which 
exhibit these features are often located outside of the clus- 
ter core, or in regions of lower X-ray luminosity. Moreover, 
the redshift evolution of the brightest cluster galaxies (BCGs) 
is peculiar and exhibits a large scatter at increasingly high 




redshift, so that at least in some cases merging could be re- 
quir ed to make them evolve into local BCGs (see for instance 
Ellis & .Tones! ll2004h ). Even if it is difficult to trace a com- 
mon "BCG evolutionary path", due to the intrinsically peculiar 
nature of th ese galaxies, some high redshift cl usters (includ- 
ing CI 1252 jRosati|l2004t iBlakeslee etai1l2003l) and C10848 
(Ivan Dokkum et alJl200ll) ) show signs of interactions or clear 
ongoing merging between few massive galaxies which could 
lead to the formation of a cD. 

However, even if the bulk of the stars had similar ages in 
the two formation scenarios (i.e. star formation occurring at 
the same early epoch in both), the epoch of assembly of the 
final mass observed locally in massive galaxies is different in 
the two cases. If merging is a relevant process in the look-back 
time range that we can probe with our observations, looking 
at progressively higher redshift one should see the number of 
massive objects decreasing as they break up into their progen- 
itors; this would cause an evolution in the shape of the mass 
function of galaxies. 

The differences in the predictions of the two formation sce- 
narios have recently become less extreme, partly due to the 
higher redshift peak of the merging activity in ACDM mod- 
els as compared to standard CDM initially considered, but also 
to different ad-hoc recipes for the star formation adopted in the 
hierarchical merging models. 

It has recently become more evident that, both 
in clu sters (references mentio ne d above) and in the 
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sive galaxies is already in place at z — 1 -H 2. However, see also 
Ivan Dokkuml (2005), finding that a considerable fraction of a 
nearby bulge-dominated galaxy sample, recently experienced 
a merging episode involving more than 20% of its final mass. 
The stellar mass function of bright galaxies shows only a mild 
evolution up to redshift 1 , cl ose to the predictio n of simple 
pure luminosity evolution (e.g. lFontana et all2004l) . 

The comparison with recent semi-analytical models how- 
ever shows that different renditions predict very different evo- 
lution, especially at higher redshift (i.e. results are very sen- 
sitive to the chosen model ingredients), and most of them 
under-produce very massive galaxies (more severely the higher 
the redshift) even when reproducing the stellar mass fun ction 
aroun d M» - however, see a lso recent results from | Bower et alJ 
(2005). At the same time, Naga mine et alJ ( 120051) show that 
with recent hydrodynamical simulations they can account for 
^ 70% of the total stellar mass at z — already being formed 
by z = 1. 

Since a direct measure of the mass function is too difficult 
at high redshift for a reasonably large galaxy sample includ- 
ing faint objects, the near-infrared galaxy luminosity function 
(LF) represents an useful cheaper surrogate. The galaxy LF is a 
first order description of a galaxy population (density of galax- 
ies as a function of their luminosity). Despite (or because of) 
its conceptual simplicity the LF has been for many years one 
of the most popular tools for the interpretation of galaxy obser- 
vations at all redshifts and in very different environments. The 
comparison of the LF at different redshifts constrains models of 
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galaxy formation and evolution ( Kauffma nn & CharloJl9 98b). 
while the comparison of the LF in low and high density envi- 
ronments probes the relevance of the environmental effects on 
the galaxy populations. The LF is historically best studied in 
rich clusters of galaxies, which provide large numbers of galax- 
ies at the same distance and, at low redshift, with high contrast 
against the background, allowing an efficient identification of 
cluster members with small contamination from background 
galaxies. At higher redshift, the faint luminosities and the sub- 
stantial background contamination makes the LF determination 
more uncertain. However, the steadily increasing data quality, 
and the quest for strict constraints on galaxy evolutionary mod- 
els, have made the study of the LFs in high redshift clusters and 
fields a popular topic. 

In this work, we determine the LF of distant (z > 1) cluster 
galaxies in the near-infrared (NIR). NIR galaxy samples are 
particularly well suited for studying galaxy evolution. Apart 
from advantages such as the smaller effect of dust extinction 
(as compared to bluer wavelengths), and the k-corrections rel- 
atively insensitive to galaxy type, they provide a relatively 
good estimate of the stellar mass in galaxies up to redshift 
z ~ 2. Therefore, near-infrared luminosity functions can trace 
the stellar mass function more effectively than bluer band LFs, 
which are more sensitive to the star formation histories of the 
galaxies. 

While LFs for cluster galaxies at low redshift (z < 0.2 + 0.3) 
have been determined for a large number of clusters, allowing 
detailed discussion of the features and the separate contribu- 
tions of different galaxy populations down to very faint mag- 
nitudes, the determination of the LF with comparable accuracy 
at high redshift is clearly more difficult. 

The NIR LF of cluster galaxies at high redshift 
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J 1998b. be Prooris et alJ 
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iKodama & Bowerl d2003l). 


Toft et alJ d2003h. lEllis & Jones 


(2004). and Toftet alJj2004). The evolution of the characteris- 



from low redshift up to z — 0.9, finding that it is consistent 
with pure luminosity evolution of a stellar population formed 
at z > 2; this result has been confirmed by subsequent studies. 
The evolution o f the faint-end slope was only studied by 
iToft et alJ d2003h and iToft et al.l l2004l> . who found a flatter 
slope at higher redshift compared to the local value. 

The adopted cosmology in this paper is Ho = 70 km 
s _1 Mpc~', Q. M = 0.3, Qa = 0.7 unless otherwise stated. 
Magnitudes are in the AB system. 

2. Data 

This work is mainly based on near-infrared images of 
three distant galaxy cl usters: RDCS J 12 52.9-2927 at z = 
1.24 (hereafter C11252. lRosati et alJ (120041. RX J08 4 8+445 3 
at z = 1.27 (hereafter C10848. IStanford et alJ dl997ll 1. 
and RDCSJ09 10+54 22 at z = 1.11 (hereafter C10910, 
Stanf ord et all l2002l) 1. The main properties of these three clus- 
ters are listed in table [l] As the data used in this work have 
already been published, we refer the reader to the papers listed 
in table|2]for details. 



We used a K, band image of the CI 1252 field obt ained 
with the ISAAC infrared imager at the VLT iLidman et alJ 
2004), a K s band image of the C10910 field obtained with 
the Pri me-Focus Infrared C amera at the Palomar 5m tele- 
scope (IStanford etal J E002). and a F160W (^ H band) im- 
age of the C10848 field obtained with the NICMOS Camera 
3 on the Hubble Space Telescope. While the quality of the 
two images for C11252 and C10848 is excellent (the PSF has 
FWHM =s 0745 and 0722 respectively, with limiting AB mag- 
nitude « 25), the C10910 image has relatively poorer quality 
(FWHM 079). While for both C11252 and C10910 the im- 
age effectively used has a radius =s 65" (i.e. slightly more than 
500 kpc in linear scale), the NICMOS image for C10848 is rel- 
atively small (the maximum radius of the mosaic is 55", i.e. 
- 450 kpc at z=1.27). For all the im ages a catalog was pro - 
duced with the SExtractor software (Berti n & Arnoutsi n"996). 
and MAG_AUTO was used as a measure of the total magni- 
tude. 

As discussed later, the galaxy luminosity function for all 
three clusters was determined by means of statistical subtrac- 
tion of the fore- and background (hereafter background) con- 
tribution. Since the images are too small to estimate the local 
background from the images themselves, a control field was 
selected for each of the cluster fields in order to determine 
the background contribution to the galaxy counts. Ideally, the 
control field should be observed in the same filter and in very 
similar conditions and depths. For C11252 the control field has 
been chosen in the FIRES (Faint Infrared Extragalactic Survey, 
iFranx et alJ feOOCh . lLabbe et all 120031) ') field in the HDF-S re- 
gion, imaged with the same instrument and in the same filter 
as the cluster field. Because of its small area, we have com- 
plemented this field at bright magnitudes with a field in the 
GOODS-S region, also observed with VLT/ISAAC (Vandame 
et al, in preparation). For C10848 the reference field h as bee n 
taken in the Hubble Ultra Deep Field (Thomps on et alJ Eo05). 
also imaged with the same instrument and in the same filter as 
the C10848 field. For C10910 we had no control field available 
imaged with the same instrument, we then selected the refer- 
ence field in the GOODS-S Ks images observed with both 
NTT/SOFI and VLT/ISAAC. Due to the similarity of the K s 
band filters used for the GOODS and C10910 images, we ex- 
pect the background estimate to be appropriate. In particular, 
the SOFI image has comparable seeing (^ 079) and compara- 
ble completeness magnitude (see figure^, and has the further 
advantage of being wider (smaller Poissonian errors). 

For the purpose of identifying point-like sources we made 
use o f the HST/ACS i mages available for a ll the cluster fields 
(Bla keslee et alJl2003l IPostman et alJl2005i) and f or the refer- 
ence regions in GOODS-S iGiavalisco et all200 4) and HUDF 
(Beckwith et al., in preparation), and of the HSTAVFPC2 for 
the reference region in the HDF-S field ( Williams^LalBoOO). 
In all catalogs, whenever possible the point-like sources were 
removed based on the MAG.AUTO vs. FLUX_RADIUS plot 
derived from the ACS images. The removed sources have 
FWHM close to the PSF of the image. Point-like sources in 
regions not covered by the ACS image were identified in the 
SOFI image itself. Since point-like sources are a small frac- 
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Table 1. Main properties of the cluster sample. Data are from Ettori et al.l J2004I) . 



Cluster 


z 


T 


R500 




M,„, 






keV 


kpc 


I0 u ergs- 1 


10 14 M o 


RDCS J0910+5422 


1.106 


6.6% 


818 ± 150 


2.83 ±0.35 


4.91 ±2.93 


RDCS J1252.9-2927 


1.237 


5-2^ 


532 ± 40 


5.99+1.10 


1.59 ±0.35 


RX J0848+4453 


1.273 


9 q+0.8 


499 ± 115 


1.04 ±0.73 


1.37 ±0.98 



Table 2. Summary of the principal characteristics of the NIR images used for the determination of the LF. Area e // is the area of 
the region actually used in this work. 



Cluster 


Telescope/Instrument 


Filter 


resolution 


completeness 
AB mag 


Area e ff 

arcmin 2 / Mpc 2 


references 


RDCS J0910+5422 


Palomar 5m / PFIC 


K s 


a/ 9 


21.5 


4.35/1.06 


Stanford et al. (2002) 


RDCS J1252.9-2927 


VLT / ISAAC 


K s 


0745 


24.5 


3.69 / 0.93 


Rosati et al. (2004) 
Lidman et al. (2004) 
Demarco et al. (in prep.) 


RX J0848+4453 


HST / NICMOS 


F160W 


0722 


25 


1.84/0.47 


Stanford et al. (1997) 

Holden et al. (2004) 

van Dokkum & Stanford (2003) 



tion of the total counts, uncertainties in their removal have little 
effect on our results. 

In order to estimate the luminosity function down to the 
faintest magnitude allowed by the data, the reference field has 
to be complete at least down to the completeness magnitude 
of the cluster field. In figure [2 we plot the number counts in 
the cluster and reference fields (normalized to the cluster field 
area). In all cases the completeness magnitude of the reference 
field is fainter or similar to that of the cluster field, therefore 
the following analysis is based on the cluster and reference 
fields catalogs down to the cluster field completeness magni- 
tude, without completeness corrections. From the turn-over of 
the number counts for objects with S /N ratio > 5, the com- 
pleteness magnitude of the cluster fields is estimated to be K s 
= 24.5 for C11252, K s = 21.5 for C10910, and F160W = 25 for 
C10848. 

This work also makes use of the extensive spectroscopic 
campaigns in thes e three clusters. We refer to [Stanford et al. 
Jl997l2002l) . lvan Dokkum & Stanford! d2003h and Demarco et 
al. (in preparation) for details on the spectroscopic follow-up 
observations. 

3. Luminosity functions 

For each of the clusters the luminosity function (LF) was calcu- 
lated by means of statistical subtraction, i.e. using a reference 
field to remove the background contribution in the assumption 
that the field galaxy density is constant all over the sky. The 
statistical subtraction of the field galaxies is often considered 
to be an uncertain method of background removal especially at 
high redshifts, where the signal of the cluster against the back- 
ground is progressively lower. However, while obtaining spec- 
troscopic redshifts for all cluster galaxies down to a reasonably 
faint magnitude is clearly unfeasible, even the determination 
of photometric redshifts relies on the availability of deep pho- 



tometry in several passbands, and on the assumption that the 
photometric redshift quality remains the same for spectral en- 
ergy distributions (SED) for which no spectroscopic redshift 
can be measured. 

Only one (CI 1252) of the three clusters studied has such 
a deep and wide photometric coverage (in addition to 38 spec- 
troscopically confirmed members), that cluster membership de- 
termination fully based on photometric redshi fts is feasible; the 
LF for CI 1252 was determined in this way bv lToft et all d2004 
(hereafter T04). Since for the other two clusters statistical sub- 
traction is at present the only viable approach, the LF of CI 1252 
was re-determined with this method, in the same region as in 
T04, as a first-order validation of the statistical subtraction pro- 
cedure in this redshift range. 

3.1. LF determination 

For each cluster the LF was determined as follows. The galaxy 
counts in both the cluster and reference fields were binned, 
and the background contribution was estimated in each bin as 
the reference field counts normalized to the cluster area. For 
each bin, the error is estimated as the sum in qua drature of the 
Poissonian errors 1 + (N + 0.75) 1/2 (Gehrels 1986) on both clus- 
ter and field counts in that bin. 

The error on the background counts should also take into 
account the effects of galaxy clustering and of the lensing mag- 
nification of galaxies beyond the cluster. However, an estimate 
of the galaxy clustering contributio n to the number cou nts er- 
ror according to the prescriptions of Huang ^etall dl997l) yields 
a negligible difference with respect to the simple Poissonian er- 
ror. Since a significant fraction of the observed galaxies are in 
the foreground of the cluster the effect due to lensing is likely 
small. Due to the large Poissonian errors we have, we can ne- 
glect these effects. 
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Fig. 1. Number counts of extended objects with S/N > 5 in the 
clusters and reference fields. Upper panel: the solid line shows 
the number counts in the C11252 field 3.7 square arcmin). 
The circles and the dotted line show the counts in the refer- 
ence regions selected in the FIRES and GOODS-S fields re- 
spectively, normalized to the cluster area. The counts in the 
GOODS-S field are shown down to the completeness magni- 
tude (K s 23). Middle panel: the solid line shows the number 
counts in the C10910 field 4.4 square arcmin). The circles 
and the dotted line show the expected background level from 
reference regions selected in the SOFI and ISAAC GOODS-S 
fields, respectively. Lower panel: the solid line shows the num- 
ber counts in the C10848 field (^1.8 square arcmin), the circles 
show the expected background level based on the Hubble UDF 
In all panels, the symbols/lines showing background number 
counts have been shifted by +0.05 mag along the x-axis to 
avoid overlapping errorbars. 



The existing spectroscopic data were taken into account, so 
that in each bin the background contribution contained at least 
as many galaxies as the spectroscopic interlopers, and the back- 
ground corrected counts were at least equal to the number of 
confirmed members. This allows a more secure determination 
of the LF in the bright end, where due to low counts (both in the 
cluster and in the control fields) the statistical subtraction may 
be poor. In the area selected for the LF determination, most of 
the bright galaxies have measured redshifts: almost 80% of the 
galaxies down to Ks = 21 for CI 1252, more than 80% down to 
Ks = 21 for C10910, and more than 85% down to F160W = 22 
for C10848. 

Since the FIRES field is quite small, in order to achieve a 
better background evaluation (and smaller errors) in the bright 
end, the background estimate for CI 1252 from the FIRES field 
was supplemented with the estimate from the control field in 
GOODS-S (ISAAC) for magnitudes brighter than 21.5. 

Due to the lower quality of the C10910 K s image (seeing 
- Of. '9), special care was taken for blended objects, particu- 
larly in the overdense cluster environment. Thanks to the avail- 
ability of HST/ACS images in passbands F775W and F850LP, 
it was possible to crosscheck the catalogs to identify obvious 
blendings. Eight cases of evident blending were identified: for 
six of them, a more 'aggressive' SExtractor configuration al- 
lowed the blended source to be split in sources located as in the 
ACS images. For the remaining 2 cases no deblending could 
be achieved, and as a zero-order approximation the flux from 
the source was split according to the flux ratio of the blended 
sources in the HST/ACS F850LP image. 

The luminosity functions are shown in figure[2] The binned 
LFs were fit with the usual [.Schechterj Jl^^ • functi on with a 
maximum likelihood method using the Cash il 19791) statistics 
C = -2£,[«,- ln(m,) - to,- ln(n,!)], where n, and to, are the 
observed number of galaxies in the i th magnitude bin, and the 
number of galaxies predicted in the same bin by a Schechter 
function of parameters M* and a, respectively. The best-fit M* 
and a are the parameters that minimize C. <t>* was not taken as 
a free parameter, but was calculated for each choice of M* and 
a by requiring that the total number of predicted galaxies equal 
the number of those actually observed. 

Even though the faint end slope cannot be well constrained 
(or is completely unconstrained, as in the case of C10910), 
due to the well known correlation of the Schechter parameters, 
leaving both M* and a free allows a better evaluation of the er- 
rors on M*. The best-fit Schechter functions are overplotted on 
the LFs in figure |2] The best fitting Schechter parameters are 
listed in table[3] 

The maximum likelihood approach gives in principle also 
an estimate of the confidence levels on the best-fit parameters, 
as if the Cash statistics is defined as above, AC is distributed 
like Ax 2 , thus AC = 2.3,6.17 gives the 1, 2-cr confidence lev- 
els for two interesting parameters M* and a. However, it should 
be noted that the Cash statistics should be applied to data which 
include background, because the background subtracted data 
are not Poisson-distributed, while the Cash statistics assumes 
Poisson probabilities. Even if we believe that the relevance of 
the 80% spectroscopic completeness on the LF bright end is 
important enough to adopt the previously described approach 
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on binned, background subtracted counts, we note that in our 
LF fitting approach we hid the fact that the counts in each bin 
have errors larger than Poissonian due to the previous statisti- 
cal background subtraction, and therefore we tend to underes- 
timate the errors on M* and a. For this reason, we also adopt a 
maximum likelihood appro ach on unbinned, no t background 
subtracted data. Recently lAndreon et al.l d2005l) summarized 
the principal reasons why one should adopt this approach, and 
proposed a method to be applied when the individual member- 
ship of the galaxies is unknown. When applying this method 
we are thus neglecting our redshift information, which means 
that we will derive conservative confidence intervals. 

For e ach of the three cluste rs we applied the method de- 
scribed in lAndreon et alJ J2005I) . taking as background dataset 
both the control regions in GOODS(ISAAC) and FIRES for 
C11252, and the control regions in GOODS(SOFI) and UDF 
for C10910 and C10848 respectively. In brief, we assume that 
the background number counts can be described by a power 
law (we used three parameters), and that the cluster LF is a 
Schechter function (also three parameters: M*, a, and <£*). We 
then find the parameters that at the same time maximize the 
likelihood for the number counts in the control field (only de- 
scribed by the power law), and in the cluster field (described 
by the power law plus the Schechter function). In all cases, we 
first searched the complete 6 parameter space for the global 
maximum, and then found the maximum likelihood on a grid 
in the M*-a plane (i.e. varying only the remaining four pa- 
rameters), so that we can draw the confidence levels for these 
two parameters. Due to very low counts (as discussed below), 
it was difficult to maximize the likelihood against M* and a for 
C10848. Since the constraints on the Schechter parameters as 
determined with the first approach are already very loose, for 
this cluster we quote them in the following. 

In figure |2 we show for reference the 1- and 2-cr confi- 
dence levels obtained with the two different approaches for 
CI 1252 and C10910. As we mentioned, the smaller ones are 
understimated but the larger ones are likely overestimated, thus 
the 'real' confidence levels are expected to lie between the two. 
The errors we quote in the following for CI 1252 and C10910 are 
derived from the larger ones. 

The CI 1252 LF is shown together with the LF based on 
photometric redshifts from T04: the two determinations agree 
within the errors, with a larger discrepancy for the last mag- 
nitude bin. The difference between the two LFs could be con- 
sidered a measure of the systematics of the two methods, with 
the error budget still dominated by low number statistics. The 
slope a determined via statistical subtraction (-0.9 + 0.3) is 
somewhat steeper than the one determined via photometric red- 
shifts (a = -0.64^ 2g), however the two estimates are consis- 
tent within the errors. The K* value (20. 1+] •£) is also found in 
agreement with the determination by T04 (K* = 20.41^? ^). 

The C10910 Schechter function plotted in figure |2] is the 
formal best fit with K* = 20.4 and a = -0.85. As it is clear 
from figure |2j the slope a is unconstrained. Assuming -1.4 < 
a < 0.4 yields an error of +1 mag on K*. 

Finally, some caveats apply to the determination of the 
C10848 LF. Namely, i) the very small field of view (less than 
200 objects brighter than i6ow — 25), resulting in signifi- 
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Fig. 2. Individual cluster luminosity functions for CI 1252 (up- 
per panel), C10910 (middle), and C10848 (lower). In all panels, 
the solid line shows the best-fit Schechter function as derived 
from the binned counts, and the inserted plot shows the 1, 2-cr 
confidence levels on the parameters M* and a, calculated from 
the maximum likelihood on binned counts (dashed lines) and 
from the maximum likelihood on unbinned, not background 
subtracted counts (solid thick lines - not shown for C10848, 
see text). The thick and thin crosses show the position of the 
best-fit parameters derived with the former and latter meth- 
ods, respectively (see text for details). In the upper panel the 
empty circles show the CI 1252 LF as determined by T04, and 
the dotted histogram shows for comparison the LF from this 
work binned with the same bin size as in T04. 
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cant Poissonian errors, ii) the likely intrinsically lower richness 
of C10848 compared to C10910 and C11252 (since C10848 has 
lower mass), iii) the presence of a known underlying superclus- 
ter in the Lynx field, and o f a lower redshift cluster projected 
in the supercluster region JStanford et alJl200li iNakata et alJ 
120051) . Mainly due to the first two reasons, the constrains on 
the Schechter parameters are quite loose, despite the signifi- 
cant depth of the F160W image. The formal best-fit Schechter 
parameters are F160W* = 20.8+} 6 and a = -1.15^. 



3.1 .1 . Composite luminosity function 

Summing up the galaxy number counts of different clusters al- 
lows the background subtraction to be more effective (averag- 
ing over uncertainties in the statistical subtraction in each sin- 
gle cluster), and the shot noise to be reduced. The composite 
luminosity function was calculated in the observed K, band at 
Z = 1.2. For this reason the F160W magnitudes at z=1.27 were 
k-corrected to K, magnitudes at z = 1.2. 

A single k-correction of 0.5 was applied at all m agnitudes, 
as derived from synthetic SEDs llBruzual & CharlotH2003l) of 
evolved simple stellar populations at that redshift (3 < age < 5 
Gyr) No correction was made for the negligible k-correction 
between z = 1.2 and z = 1.106 or z = 1.237. All the indi- 
vidual LFs were binned with a bin size of 0.5 mag (binning 
was adjusted to optimize individual magnitude coverage taking 
into account completeness limits and corrections for different 
redshift). C11252 and C10848 have the same "K," band com- 
pleteness after these corrections, so their LFs were just summed 
up, and the errors were added in quadrature. C10910 instead is 
much shallower, so its LF was added to the composite LF up 
to its completeness magnitude. The composite LF beyond this 
magnitude is computed as the composite LF without C10910 
multiplied by the ratio of the total counts (including all three 
clusters) to the total counts excluding C10910, computed in the 
magnitude interval where the C10910 photometry is complete 
(and errors were scaled accordingly). Due to the bright com- 
pleteness limit for C10910, and the low counts for C10848, it 
is clear that the faint end of the composite LF is dominated 
by C11252. In our case, building the com posite LF follo w- 
in g other comm on methods as described in lcollessl (fl989) or 
in G arilli etafl dl999h . produces results consistent within 1-cr 
w ith the LF calculated as above (note that the method described 
in G arilli etafl l ll999 ) tends to give smaller errors, and a flatter 
slope than the one in [collessl(ll989h V 

The composite LF was also derived with a maximum like- 
lihood approac h on unbinned, not bac kground subtracted data, 
as described in lAndreon et al.l l|2(X)5). However, since the LF 
of C10848 is measured in a different passband than those of 
C11252 and C10910, C10848 was not included in this case. The 
composite LF was determined by fitting at the same time the 
counts in the two cluster fields (CI 1252 and C10910) and in all 
the Kj band control fields. The two cluster LFs are assumed 
to be described by the same M* and a (but clearly have two 
different <£*). 

The composite luminosity function is shown in figures [3] 
and |H The best-fit Schechter has K* = 20.5±° 4 and a = 
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Fig. 3. The composite cluster luminosity function at z=1.2 
(filled dots) compared to the field galaxies LF at z - 1. The 
solid line shows the best-fit Schechter, and the inserted plot 
shows the 1, 2-crconfidence levels on the Schechter parame- 
ters as in figure |2 All field galaxies LFs have been arbitrarily 
rescaled. 
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Fig. 4. The composite cluster luminosity function at z=1.2 
(filled dots) compared to the local cluster galaxies LF, corrected 
by 1.3 mag for passive evolution (see text). The solid line shows 
the best-fit Schechter, and the inserted plot shows the 1, 2- 
crconfidence levels on the Schechter parameters as in figure|2| 
All local cluster galaxies LFs have been arbitrarily rescaled. 



-1.0 + Q 3. In fi g ure [3] w e also show t he fiel d LF at z ^ 1 from 



IPozzetti et all J2003h . iDrorv et alJ J2003) and bahlen et ail 
(2005). In figure |4] the composite LF is compared to the lo- 
cal cluster galaxies K band LF (corrected by 1.3 mag for 
passive evolu tion as deriv ed b elow) as measured in C oma 
bv lDePropris et alJ Jl998l) and lAndreon & Pellol lE)00), and 
in sampl e s of nearby clusters by Balo gh et alJ (J2001) and 
iLin et alJ ([2004). All the local and field LFs have been arbi- 
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Fig. 5. The background number counts as estimated in differ- 
ent control fields. The filled symbols show the counts from the 
control field in the FIRES region, the solid and dashed lines 
from those in the GOODS-S regions (ISAAC and SOFI, re- 
spectively), and the dotted lines show the counts from the 17 
ISAAC tiles in the GOODS-S field. 



trarily rescaled. Both the Coma LFs shown were measured in 
the H band and shifted to K band w ith a colour term H - K - 
0.24. The K band field LFs from IPozzetti et all fe003l) and 
IProrv etalJ j2003), and the J band field LF from bahlen et ail 
J2005h . were converted to observed K, ma gnitudes at z = 1 .2 
by equations 2 and 1 in P ozzetti et alJ J2003b . 

With these data, no significant difference can be seen be- 
tween the shapes of the cluster and field luminosity functions 
in the probed magnitude range, even if we find some evidence 
of an excess of very bright galaxi es with respect to the field, as 
suggested for instance in lDe Propris et all d2003l) . 

The derivation of luminosity functions based on statistical 
subtraction may be affected by field-to-field variations, as the 
background evaluated from a control field may not be represen- 
tative of the background in the cluster field. A robust estimate 
of such an effect requires adequately deep and large K band 
fields. The VLT/ISAAC observations of the GOODS-S field, 
covering ~ lOOarcmin 2 , are currently the best data set avail- 
able for this purpose. We note that part of this same field has 
been used as control field for determining the LFs of C10910 
and C11252. In each of 17 selected ISAAC tiles in GOODS-S, 
we considered all galaxies within the central region (with the 
same area as the CI 1252 field used for the LF), down to the 
photometric completeness. We also considered in this case as 
an "additional tile" the FIRES HDF-S field we already used as 
control field for C11252. 

In figure|5j we show the background number counts as esti- 
mated from different control fields: those from the FIRES field, 
which was used as control field for CI 1252 (area 5 arcmin 2 ) 
at faint magnitudes, those from the ISAAC GOODS-S mosaic 
used to complement the C11252 control field at bright mag- 
nitudes (area 53 arcmin 2 ), those from the SOFI GOODS-S 
mosaic used for C10910 (area ^ 152 arcmin 2 ), and those from 
the 17 ISAAC tiles in GOODS-S (area * 3.7 arcmin 2 each). 



We then built the C10910 and CI 1252 LFs as described 
above (making use of the spectroscopic information), based on 
these 18 small control fields, and for each of these we redeter- 
mined the composite LF. For C10848 we always used the LF 
determined above, since due to its very large uncertainties it 
has lower weight compared to C10910 and C11252. Thus, we 
obtained 18 composite LFs corresponding to different control 
fields. 

The GOODS-S ISAAC data, while being considerably 
wide, are not deep enough to reach our faint-end magnitudes. 
Therefore, we complemented the counts from each ISAAC tile 
at magnitudes fainter than its completeness with the FIRES 
counts. In this way, we probe field-to-field variations at magni- 
tudes brighter than K s ^23. The medians and standard devi- 
ations of M* and a obtained from these 18 determinations are 
M* = 20.5 mag, cr M .=0.14 mag, a = -1.0, cr a =0.06. 

In order to also account for field-to-field variations at 
fainter magnitudes (probed only by FIRES), we repeated the 
LFs determination by normalizing the FIRES counts at faint 
magnitudes by the ratio of the ISAAC/FIRES counts at magni- 
tudes brighter than the completeness in each tile, assuming that 
the number density ratio is the same at fainter magnitudes. The 
medians and standard deviations of M* and a obtained from 
the 18 redeterminations with this procedure are M* = 20.5 mag, 
o- M .=0.09 mag, a = -1.0, cr a =0.09. 

In figure[6] we show the results from these two sets of tests. 
In the main panel, we show the confidence levels as shown in 
figures|3]and|4] and we overplot the Schechter parameters M* 
and a obtained for the 18 LFs. Since most of the data points 
overlap near the original determination marked by the cross, 
we show their distributions in the side-panels. 

3.2. Evolution of the restframe K band LF 

It is customary to compare the observed evolutio n of M* with 
redsh ift with different predictions. Even if initially Bar ger etafl 
dl998h found n o significant evolution f or the infrared M* up 
to z > 0.5, first D e Propris et alJ ( 1999) an d then other works 
on hi 
| 2003 [ 

2004) found that the evolution of M* up to z ss 1 is con- 
sistent with pure luminosity evolution of the cluster galaxies, 
an d inconsi s tent w ith no-evolution pr edictions. As discussed 
in lAndreonl d200l|) and lAndreonl J2004I) . measuring an evolu- 
tion in the LF from a change in the best-fit Schechter param- 
eters is not straightforward. A density (clustercentric radius) 
dependent LF would imply a dependence of M* (and a) on 
the surveyed area in different clusters at different redshifts, 
and the correlation of M* and a could introduce spurious re- 
sults. This makes it difficult to study the LF evolution exclu- 
sively based on the evolution of the characteristic magnitude 
M*. However, it is unlikely that the bright end of the LF is 
dominated by galaxies in the cluster outskirts, and even though 
the LF has indeed been found to be dependent on the sam- 
pled region within the cluster, this dependence mostly affects 
the LF at magnitudes fa inter than those we can probe in this 
work (e.g.. lPopesso et af] J2005)). When comparing LFs from 
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Fig. 6. Impact of field-to-field variations on the determination 
of the LF. In the main panel the best-fit M* and a (cross) and 
their 1- and 2-cr confidence levels are shown, as in figures [3] 
and|4] The M* and a parameters from LF determinations with 
18 different small control fields are shown as squares. In the 
side-panels the distributions of these M* (top) and a (right) 
are shown. Solid symbols (in the main panel) and solid lines 
(in the side panels) show the results for GOODS-S tiles counts 
complemented with FIRES counts at faint magnitudes, while 
empty symbols and dotted lines show the results with modified 
FIRES counts at faint magnitudes (see text for details). 



different studies, we note that our LF is based on the central 
r < 500 kpc clust er regions, approxi mately corresponding to 
r^oo, while the De Propr is et alJ i 19981) LF is determined in the 
central r < 350 kpc Coma region, the lAndreon & Pellol feOOOh 
LF is determined in a ^ 500 x 500 kpc 2 regio n offse t by 360 
kpc from the Coma centre, and the lLin et alJ (12004) LF is de- 
termined within the virial radius and is found to be very similar 
to the LF derived within t^qq. 

Finally, if a is free in the Schechter fit, the errors on M* are 
reliable enough to make a fair comparison of the M* evolution 
with different predictions. 

The comparison with previous determinations of K* at 
lower and similar redshift is shown in figure [7] In agreement 
with previous work, the measured K* is consistent with passive 
evolution predictions for an L* galaxy formed at z> 2. 

Converting the observed K* to the absolute K, band magni- 
tude via K, t = K s „ hx - 51og(d L /10 pc) - (K S)rest - K s , obs ) z , as in 
IPozzetti et all fe003l) . gives Kt r „„ = -23.41 + ° 4 . C ompared to 
the Coma LF K* ^ -22.15 (De Propri s et al J 19981) . this yields 
an evolution of AK* = -1.3*j' . As shown in figure^ the shape 
of the composite LF is very similar to the shape of the local 
cluster galaxies LF shifted 1.3 magnitudes brighter. 

This shape may be quantified in a non-parametric way 
by an analogue of the 'giant-to-dwarf ratio' (GDR), which 
is defined in this case as the ratio of the number of galaxies 
brighter than K s = 21.2 to the number of galaxies with 21.2 < 
Kj < 24.2. The K s = 21.2 threshold corresponds to an abso- 
lute magnitude of -22.7, and to a stellar mass of 8 x 10 10 
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Fig. 7. The redshift evolution of the characteristic magnitude 
K*. Different symbols show different determinations of K* as 
indicated in the legend. The no evolution prediction is calcu- 
lated from the Coma K* toe Propris et al.ll 19981) , k-corrected 
as in IPozzetti et alJ J2003I). Passive evolu tion model predic- 
tions are from Kodama & Arimoto normalized to the 
Coma K* . All d e terminations have fixed slope a = -0 .9, except 
lAndreonl feOOlh . iToft et alJ fc003fc. iLin et alJ J2004> . T04 and 
the present work. For lLin et al.l (2004) the plotted value is the 
mean of two values determi ned for a = -0.85 and a = —1.1. 
For both this work and lToft et alJ J2003I) also the K* value with 
a fixed at -0.9 is shown (corresponding empty symbol). Note 
that errors on M* which are computed with a fixed and free are 
not directly comparable. Error bars on the x-axis, when plot- 
ted, represent the redshift range of clusters which are combined 
to draw that point. 



M (for a Salpeter IMF). For our composite LF this GDR is 
0.2 + 0. 1 . If we estimate this GDR with the same absolute mag- 
nitude cut using the Coma LFs as shown in figure|4](i.e. keep- 
ing into account the 1.3 mag brightening), we similarly find a 
GDR of 0.2 - 0.3. This suggests that a large fraction of the gi- 
ant population was already in place at z ~ 1. We remind that a 
computation of evolution corrections, k-corrections, and H-K 
colour terms at redshift zero, are involved in such a compar- 
ison. We also remind that we are applying a single 1.3 mag 
brightening for the whole LF, which is clearly a simplistic as- 
sumption, since galaxies with different star formation histories 
have different evolution corrections. 

3.3. Contribution to the LF from early and late type 
galaxies 

The study of the early and late type galaxies LFs at this redshift 
is challenging. 

Even when deep ACS data are available,(as in the case of 
CI 1252), morphological analysis is not feasible at the faint end 
(zab <: 25 ). Here we use the m orphological information pre- 
sented by (Blakeslee et al. 2003), which is based on i or z band 
ACS imaging. Due to the red B - z restframe colours (corre- 
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sponding approximately to the observed z-K), it is not possi- 
ble to build a morphological catalog down to the K, band com- 
pleteness magnitude. On the other hand, if early and late types 
are distinguished based on their colours, while the cluster red 
galaxies can still be isolated with relatively small background 
pollution, this does not hold for the blue population. For these 
reasons we only attempt a separation of the contributions to the 
LF from early and late type galaxies for the bright end of the 
C11252LF. 

For this cluster we can use the photometric members se- 
lected in T04, and distinguish early and late types based on 
the best-fitting template from photometric redshifts. The pho- 
tometric redshifts in T04 were determined with 7 passbands 
again st four templates (E/SO, Sbc, Scd, Irr) fromlColeman et alJ 
( 1980), two starburst SEDs from lKinnev et alJ i 19961) . and in- 
terpolations between these six SEDs (see T04 for details). We 
then defined as early-types those galaxies best fitted with SEDs 
earlier than midway between E/SO and Sbc (i.e. roughly includ- 
ing E/SO and possibly some Sa galaxies). We can then sepa- 
rate the two contributions to the LF, virtually down to the K s 
band completeness magnitude, using the completeness func- 
tion calculated by T04 (their fig. 5) to correct both the early and 
late-types LFs for incompleteness due to photometric redshifts. 
While the reliability of photometric redshifts may be uncertain 
at the faint end, the early/late-type separation at the bright end 
is robust and the completeness correction negligible. 

In figure~8]the early (filled dots) and late-types (empty cir- 
cles) bright end LFs are shown, where the separation in early 
and late-types reflects their SED properties. We also show as a 
dashed li ne the LF of ear l y type galaxies morphologically se- 
lected bv lBlakeslee et all 120031) (brighter than K s < 22.5, as 
the typical z-K colours drive the sample beyond the Zg5o com- 
pleteness limit at fainter magnitudes). The morphologically and 
SED-selected early-types LFs are in very good agreement, 
suggesting that the bright end of the LF is already dominated 
by early-type galaxies, either selected on their morphology or 
on their spectrophotome trie properties. 

A histogram of the red-sequence galaxies (which are ex- 
pected to be mostly early types) is also shown for comparison 
This was derived taking all the galaxies (within 65'.' from the 
cluster center) with 1775-Z850 colours 0.16 mag redd er or 0.14 
mag b luer than the red sequence determined by Blak eslee et alJ 
(2003). This is mu ch larger than the intrinsic scatter found 
by blakeslee et alJ (120031) . however the colours we used are 
1'.'5 aperture colours, correspondingly the scatter is expected 
to be larger. No statistical subtraction was attempted in this 
case, since a reference field with deep enough K s , 1775 and Z850 
imaging is available, however spectroscopic interlopers were 
removed. As a result, the histogram shown is an upper limit 
to the effective LF of red sequence galaxies. The shaded area 
represents the 16-84 percentile variations of this histogram due 
to photometric errors, and was derived by simulating 100 cat- 
alogs where the i775-Z85o colour was randomly shifted within a 
Gaussian of cr equal to the photometric error on the 1775 - zgso 
colour. This histogram also confirms that the LF bright end is 
largely dominated by galaxies hosting evolved stellar popula- 
tions. 
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Fig. 8. The contribution of early and late type galaxies to the LF 
of C11252. The filled and empty symbols show the LFs for early 
and late type photometric members (based on the T04 sam- 
ple), classified from their broad band colours. The dashed line 
shows the contribution of morphologically selected early-types 
(Blakeslee et al. l2003h . The solid histogram shows the num- 
ber counts of all galaxies (excluding spe ctroscopic in terlopers) 
along the red sequence as determined in (Blakes lee et al.l2 003) 
(no statistical subtraction of the field galaxies contamination 
was made; the shaded area shows the effect of photometric er- 
rors - see text). The solid curve shows the best-fit Schechter 
function for the CI 1252 global LF. 



A solid determination of the LF of red-sequence galaxies 
would require even more extensive redshift information at the 
faint end, which is however beyond the current spectroscopic 
limit. With our data, we observe some evidence of a deficit of 
faint galaxies on t he red-sequence, wh i ch has been reporte d 
in other studie s JKaiisawa et all J200QI) . iNakata etafl J200ll) 
iKodama et aD J2004I) . be Lucia et alJ J2004I) . iTanakaetalJ 
il2005l) 1. and is usually interpreted as a sign of downsizing. 



3.3.1 . Contribution from clusters to the bright 
galaxies budget 

A rough estimate of the contribution of the bright cluster galax- 
ies (< M* +2) to the total bright (< M* +2) galaxy budget can 
be obtained, by combining the cluster galaxies LF (and its mea- 
sured evolution) with the known space density of clusters out 
toz ~1. 

Our results show that at least at magnitudes brighter than 
M* +2, the cluster galaxies LF appears to evolve mainly by 
passive evolution up to z — 1 . In addition, the normalization of 
the K band galaxy LF of X-ray luminous clusters was found 
to be consistent with the loc al one at least out to z - 0.8 
jTrentham & Mobasheilll 99*81) . Therefore, we can assume that 
high-redshift clusters of a given X-ray luminosity contain sim- 
ilar numbers of bright (< M* +2) galaxies as low-redshift clus- 
ters of the same X-ray luminosity. 
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We use the relation between the number of galaxies brighter 
than M* +2 within r2oo, and the cluster mass within r2oo, as 
derived by Popesso et al. (in preparation): 



N ga i axies - 10- 712 (M 200 /M o ) - 61 



(1) 



and the relation between the X-ray lumi nosity and the clus- 
ter mass within r2oo (Reiprich & Bohringer 2002), in the form: 



log 



^X{().\-2AkeV) 



/i^l0 40 ergs- 



-19+ 1.58/og 



( M 200 ) 



V'*50 



(2) 



By combining these two scaling relations (and converting 
to Ho = 70 km s~' Mpc~\ and L x (0.5 - 2.0/teV)[10 44 erg s" 1 ]), 
one obtains the number of bright galaxies in a cluster of lumi- 
nosity L x : 



ML X ) 



74.5 • (L x (0.5 - 2.0/fceV)[10 44 erg s" 1 ]) 386 . 



(3) 



In the assumption that this holds in the redshift range [0+ 1 ] , 
we can use the redshift evolution of the cluster X- ray lumi- 
nosity function jRosati et al.ll2002l iMullis et aDl2004 to esti- 
mate the number density of clusters at a given luminosity and a 
given redshift. By using the N{L X ) relation above, we can com- 
pute the number density of bright galaxies in the cluster virial 
regi ons at each redshift. A comparison with th e field LF (we 
used lTrentham et alJ d2005l) at z - 0, and lDahlen etaP d2005h 
at z - 1) yields the contribution of cluster galaxies to the total 
bright galaxy budget. We thus find that a fraction of less than 
10% (— 6 -i- 7%) of the bright galaxies at z — is located in the 
virial regions of X-ra y luminous clusters (in broad agreement 
with what reported in lDe Propris et all d2003h '). and similar (^ 
5%)atz^l. 

4. Mass to light ratios 

The cluster dynamical mass-to-light (M/L) ratio has long 
been a matter of interest due to the fact that clusters collapse 
from regions several Mpcs wide, thus carrying both their 
mass and galaxy content from representative portions of the 
Universe. Their M/L ratio could thus be similar to that of the 
whole Universe, even thought this is much dependent on how 
different galaxy evolution is in such high density environments 
with respect to average density regions. 

The cluster M/L ratio can in principle be used to study 
how galaxies evolve in dense environments. The M/L ratio is 
known to in crease with the syst em mass going from galaxies 
to clusters (Bahc all et al.lll995l) . in agreement w ith the pre- 
dictions of models of bia sed galaxy formation (tDavis et alJ 
ll985UBardeen et al.lll986h . However, it is not yet completely 
clear, also in recent studies, whether a mass dependence 
of the M/L ratio is observed in the limited mass range of 
groups and clusters. While some works find a measurable in- 
crease of the M/ L ratio with the t ot al mass of t h e system 



Girardietal. (20Q(j), 


Hoekstra et alJ 




(2001), GirardLeLalJ 


(2002), 
(2002), 


Marinoni & Hudson 
LinetalJ d2003l 120 


d2002l). [Bahcall & Comerford 
04h. iRinesetalJ d2004k. iTullvl 



2005)), some others conclude that the M/L ratio is ap- 
proximat e ly the same in gro u ps and clusters (for instance 
Dresslerl Jl978ll. lDavidet al.1 dl995l). ICarlberg et alJ (l996 



uresslerl '"avid et al. 

1997), Ci rimeleetal] lM997). : 



SEooo)). 



As observed for instance in Tullv (2005), galaxy groups 
with lower B band M/L ratios compared to more massive sys- 
tems generally have a population of late type galaxies with 
on-going star formation. Moreover, systems with predominant 
E/SO/Sa population tend to have larger M/L ratios; this would 
be consistent with more dense regions forming earlier, and 
therefore ending their star formation earlier. More recently, the 
K band M/L ratio, having very little dependence on on-going 
star formation, has been found to be mass-dependent dLin et alJ 
2003 I2004T) . with K band light per unit mass being higher by a 
factor =s 2 in low mass clusters th an in massive ones. 

As discussed for instance in iRines et alJ (|2004), the ob- 
served mass dependence could be due either to processes like 
tidal stripping and dynamical friction disrupting galaxies in 
massive clusters, or to a reduced star formation efficiency in 
such systems (possibly the heating of the ICM cutting of f the 
supply of cold material needed to form stars dBlanton et alJ 
ll999NBaloghetall2000l) 1. 

Keeping in mind these issues that complicate the compari- 
son of inhomogeneous samples, as well as possible systematic 
differences i n the masses estimated by di ff erent means (see 
for in stance LSanderson & PonmanI J2003I) . lAndernach et alJ 
d2004l) '). we can compare the M/L ratio of the clusters stud- 
ied in this work with estimates at lower redshift. A compar- 
ison of the cluster M/L ratios in the B band, out to redshift 
z — 0.8, was presented for instance bv iHoekstra et al.l (|2002), 
who found that the evolution of the M/L ratio is consistent with 
the luminosity evolution of galaxies as derived from the funda- 
mental plane in distant clusters (see their figure 14). 

We have estimated the K band M/L ratios of C10910, 
CI 1252, C10848 making use of th e LFs derived abo ve and of 
the X-ray mass profiles derived in lEttori et al.N2004l) . 

Once the Schechter parameters a and K* est (and the corre- 
sponding characteristic luminosity L*) have been determined, 
the K band projected total luminosity within the surveyed area 
can be calculated via direct integration of the Schechter func- 
tion (L, , = <&*LT(2 + a)). The K band luminosities for the 
three clusters within the surveyed area (as listed in column 6 of 
table|2j are listed in table[3] We assume K = 5.2 when calcu- 
lating luminosities in units of L . The errors were determined 
by calculating the luminosity with Schechter parameters run- 
ning on the 1-cr confidence level for K* and a. 

The projected M/L ratios for the three clusters were derived 
within t he surveyed area, u sing the projected mass profiles de- 
rived in lEttori et all 12004) (see tableQ. 

Note that the mass to light ratio derived for CI 1252 by sum- 
ming up the luminosities of the photometric members from T04 
without further corrections (M/Lk = 17^M O /L ) is consistent 
with the quoted value. By considering only the confirmed spec- 
troscopic members, we obtain M/Lk — 31 + 5M e /L Q , which 
should be considered as an upper limit. 

It should be noted that, since our surveyed areas are differ- 
ent from each other and are small, we rely on the assumption of 
a negligible or small dependence of the M/L ratio on the clus- 
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Table 3. Luminosity function parameters and estimated absolute K band luminosities and mass-to-light ratios. Column 2: in- 
dicative radius within which the LF and the M/L ratio are measured; column 3: original passband in which the LF is measured; 
column 4: LF characteristic magnitude M* as measured in the original passband, column 5: M* k-corrected to the restframe K s 
band (F160W was previously corrected to K s band by a factor -0.5); column 6: LF faint end slope a (note that for C10910 a is 
unconstrained, and the error on M* for this cluster is estimated assuming that -1.4 < a < -0.4); column 7: the total restframe K s 
luminosity within the effective area listed in column 6 of table|2j column 8: the restframe K s band mass-to-light ratio (the errors 
come from the errors on the projected mass and on the total luminosity). Note that for C11252 we report the total luminosity and 
the M/L ratio derived from both the LFs from this work and from T04 (first and second value respectively). 
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tercentric distance within the virial radius fanes et alJl200ll 
iKneib et alJ2003llRines et alJ2004 . when comparing the M/L 
ratios derived here with other measurements. 

In figure [9] the K band M/L ratio for the three clus- 
ters is compared to previous deter mination s at lower red shifts 
JCarlberg et al.lll997llc1rimele et alJll997l Girardi et al. 2000, 



Hradeckv et al. 2000, 


Rinesetal. 2001, 


Hoekstra et al. 2002 


Kneib et al|2003[|Lin et al|2003USanderson & Ponman 2003 



Andernach et all2004lGavazzi et all2004llRines et alJ 2004). 

The M/L ratios published in passbands different from K 
were rescaled to the K band using the colours of a simple stellar 
population formed at z=5 and the AB colours of the Sun (B - 
K) =O.15, (R - K) = -0.76, (V - K) B = -0.33. These rescaled 
measurements are shown with empty symbols. 

To avoid excessive confusion in the plot, some points 
do not represent a single cluster but are based on different 
samples: 32 groups and clusters from Sanderson & Ponman 
l l2003l) . 8 group s and clusters from [ Hradeckv et al J J2000). 
16 clusters from ICarlberg et all i 19971) gro uped in 4 redshift 
bins, 105 clusters fro m lGirardi et alJ ll2000l). 12 c l usters from 
ICirimele et all ll 19971). 13 clusters fromlLin et alJ lEool . 180 
clusters from lAndernach et alJ ll2004 . and 9 clusters from 
iRines et alJ fe004 . In such cases, the weighted average (and 
corresponding error) of the sample is plotted in figure|9] 

Some of the M/L ratios plotted in figure|9]were derived with 
masses estimated from kinematics (crossed points) or from 
strong/weak lensing (circled points). The solid line traces the 
expected evolution of the M/L ratio, neglecting any evolution in 
the dark halo mass, in the assumption that the M/L ratio evolves 
following the luminosity evolution of the cluster galaxies, for a 
pure luminosity evolution of a simple stellar population formed 
at z — 5 (normalized at M/L = 51 at redshift zero). 



5. The stellar mass function 

At redshift 1, the K s band (rest-frame 1/vm) luminosity is 
still a good tracer of the stellar mass. Therefore, we can draw 
an estimate of the stellar mass function (MF) of cluster galaxies 
at redshift z =£ 1.2 from the composite K s band LF. The K s band 
light translates into the stellar mass via the stellar M/L ratio of 




0.5 1.0 
redshift 



1 .5 



Fig. 9. The evolution of the restframe K band mass-to- 
light ratio. The filled circles are the M/L ratios from this 
work. The star shows the M/L determined for CI 1252 using 
the T04 L F. The filled squar e s are K band M/L measure- 
ment fr om IRines et alJ (l200ll) . IKneib et all J2003I) . Lin et alJ 
( 20031). IRines et alJ ll2004h . The empty symbols show the K 
band M/L ratio derived from published M/L estimates in 
other p assbands, i.e. B/B , b a nd (circle) from | Girardi 1 
feOOd). iHoekstra e t al. (20ol. ISanderson & Ponmanl J2003 
iGavazzi eUriTlToO^. R/r ban d (squa re) from ICarlberg et alJ 
dl997h andt Andernac h et al] ll2004. and V band from 
ICirimele et al.1 (1 19971) and iHradeckv et al.1 J2000I) (see text). 
Crossed and circled points have masses determined from 
galaxy kinematics or strong/weak lensing respectively, at a dif- 
ference with X-ray masses used in this work. The solid line 
traces the evolution expected from pure luminosity evolution 
of a simple stellar population formed at z=5, assuming that the 
cluster M/L ratio only evolves because of stellar evolution in 
galaxies, and neglecting any evolution in the dark mass. The 
dotted line shows the reference (z=0) value. 



each galaxy's stellar population, which depends on the galaxy 
star formation history (SFH) and on its age. 

In principle, if we knew which galaxies are contributing 
to the composite LF (and we measured their photometry in a 
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sufficient number of passbands) we could measure the stellar 
M/L ratio for each of the galaxies via SED fitting, and directly 
determine the stellar mass function. However, since the com- 
posite LF was determined in a statistical fashion, we also have 
to statistically evaluate the stellar M/L ratios of cluster galax- 
ies along the LF. To this aim, we used the CI 1252 photomet- 
ric members selected by T04 to statistically estimate, in each 
magnitude bin of the LF, the median and scatter (16'''-84' /! per- 
centiles) colours iws-Kj, Z85o-K s , and J s -K s . In the assumption 
that the contributions to the LF from different galaxy popula- 
tions in C10910 and C10848 are approximately similar to those 
in CI 1252, these colours obtained along the CI 1252 LF will be 
representative of those along the composite LF. 

We then built a set of 160 synthetic SEDs with the 
iBruzual & Charlotl J2003h code, with (delayed) exponentially 
declining SFH with 0.005 < r < 1.5 Gyr and 0.3 < age < 5 
Gyr, solar metallicity, Salpeter IMF and no reddening. For each 
magnitude bin, we selected all the models with the appropri- 
ate colour (median + scatter as above) for that bin, thus de- 
termining a rough estimate of the M/L ratios in that magni- 
tude interval. For each of the three colours i775-K s , Z85o-K s , 
and J s -K s we extracted 100 realizations of the MF perturbing 
the number of galaxies in each magnitude bin with a cr equal to 
its Poissonian error and then spreading in mass this perturbed 
number of galaxies within the above defined M/L ratio range 
for that magnitude bin. We then considered the median and 
minimum/maximum MFs over the 100 realizations, for each of 
the three colours, which give three estimates of the MF which 
are perfectly consistent. 

We thus averaged these three estimates obtaining the MF 
plotted in figure [H)] the shaded region corresponds to the aver- 
age minimum/maximum MFs computed according to the above 
description. We also show for comparison the MFs for field 
galax ies at similar redshift JDrorv et al J l20o4 iFontan a et alJ 
2004), for field galaxies at redshift zero (Cole et al. 2001), and 
for local cluster galaxies (Balog h et al.1120011) 1 . all arbitrarily 
rescaled so to match at M stars 5 ■ 1O 1O M . The MF is shown 
for masses greater than 1O IO M , which is the estimated 
mass completeness determined for our completeness magni- 
tude, considering the M/L ratio of a stellar population formed 
at z^lO with subsolar metallicity and no reddening. The stellar 
mass corresponding to the characteristic magnitude K* = 20.5 
is approximately 10 11 M . 

Within the uncertainties affecting this MF determination, 
the shape of the MF of massive objects at redshift 1.2 is not 
significantly differ ent from the local one as measured from the 
2dF galaxies JCole et alJkOol . This is in agreement with the 
very mild evolution of the mass function observed in the red- 
shift range [0:1] bot h in the field (e.s 



,g., FontanasLalJ fctXMh 
verl (|2003)). The presence 



and in clusters ("e.g.. lKodama & Bower! 
of massive objects (M sta r S > 1O U M ) is independently con- 
firmed at least in CI 1252 from SED fitting on 9 passbands from 
restframe NUV to NIR (Rosati et al., Rettura et al., in prepara- 
tion). However, we note that by probing galaxy clusters we are 
only sensitive to the evolution of the shape of the MF, and not 





■, 
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cluster galaxies at z— 1.2 

local cluster galaxies (Balogh et al. 2001) 

AA field galaxies at 1<z<1.5 (Fontana et ol. 2004) 
■ field galaxies at 1<z<1.2 {Drory et al. 2004) 
O local field galaxies (Cole et al. 2001) 




The MF from B alogh et alJ l 200 1 ) was calculated with a Kennicutt 



IMF, so we applied a correction logM™^'" 



: logM* 



+0.35. 



Fig. 10. The stellar mass function as determined from the com- 
posite K, band LF. The shaded area accounts for errors in the 
K s band LF and in the stellar mass to K, light ratio (see text for 
details). Other determinations of the stellar mass function are 
shown for reference, as indicated in the legend. 



to its normalization, whereas field galaxy survey s find the evo- 
lution of the MF to beco me significant at z > 1 JFontana et alJ 
l2004lDrorv et all2004l) . 



6. Conclusions 

We have studied the near infrared luminosity function of high 
redshift cluster galaxies in three X-ray luminous clusters that 
are among the most distant discovered so far (1.1 < z < 1.3). 

These clusters bear the strongest leverage on evolutionary 
studies as they probe a redshift regime when the Universe was 
less than half of its present age. By measuring the K band lu- 
minosity function (and galaxy stellar mass function) in these 
systems, and by comparing it with that at z — 0, we can set 
valuable constraints on the galaxy evolution in dense environ- 
ments in the redshift range [Oh- 1 ] . 

We should note, however, that derived quantities still suf- 
fer from large uncertainties, due to small-number statistics and 
possibly biases in the galaxy populations (due to probing only 
the cluster cores). 

The LF resulting from this work is consistent with previous 
determinations at similar or lower redshifts, and consolidates 
a scenario where the evolution of the characteristic magnitude 
M* (AM* a; -1 .3) is consistent with predictions of passive evo- 
lution for a stellar population formed at z > 2. Moreover, we 
find that the overall shape of the high redshift LF matches the 
one of the local cluster galaxies LF, once such a brightening is 
taken into account. 

Similarly, the evolution of the K band LF of field galaxies 
has been found to be consistent with passive evolution up to z = 
1, with a density evolution lo wer than 30% and a brighteni ng of 
AM* ^ 0.5 -r 0.7 mag (e.g IPozzetti et all d2003h . iDrorv et alJ 
d2003l) - however, see also Dahle n et al.N2005h . 

A direct comparison of the field and cluster LFs at redshift 
zero, has reveal ed a significant diffe r ence both in the B and i n 
the K band (e.g. Balog h et alJ d200ll) . be Propris et alJ (E)03)1 
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amounting to 0.3 mag in M* and =2 0.1 4- 0.2 in a. While 
there is a hint of exceeding very bright galaxies in our distant 
clusters with respect to the field, our error bars do not allow us 
to make such a claim. 

For CI 1252, for which best quality data and extended wave- 
length coverage are available, we attempted a separation of the 
contributions of late and early type galaxies to the bright end 
of the LF. We find that, already at z — 1 .2, the bright end of the 
LF appears to be dominated by early type galaxies selected ei- 
ther on the basis of their morphological appearence or of their 
spectrophotometric properties. 

Using the individual LFs for the three clusters, we calcu- 
lated the K band cluster M/L ra tios m aking use of X-ray mass 
profiles derived in Ett ori et all (J2004). The M/L ratio tends to 
be smaller than the typical value at redshift zero, as expected 
on the basis of pure luminosity evolution of the cluster galax- 
ies stellar populations. However, a much larger sample would 
be needed for a detailed investigation of the evolution of the 
cluster M/L ratio. 

Finally, from the composite K, band LF we have estimated 
the stellar mass function of cluster galaxies. The observed K s 
band light at z — 1.2 corresponds approximately to the rest- 
frame z band light, and is considered as a good tracer of the 
stellar mass. We have outlined in the introduction how the de- 
termination of the stellar mass function at high redshift sets 
strong constraints on the evolution of massive galaxies. While 
the early formation epoch of the bulk of the stars in massive 
galaxies is now generally established by several observations, 
the epoch of the major mass assembly can only be assessed by 
studying the redshift evolution of the mass function. Our study 
shows that the massive (M stars > 1O 1O M ) galaxy populations 
in massive high redshift clusters have not significantly changed 
since z — 1, apart from passive evolution of their stars, thus 
extending previous results at lower redshifts. The shape of the 
stellar mass function at z — 1 .2 is found to be consistent with 
the one observed in local clusters, within our uncertainties. 

The high-mass end of the LF, made of giant (M stars > 
1O H M ) E/SO galaxies, is already in place at z — 1.2. This 
points toward an early assembly of the galaxy mass, mostly 
completed before z - I, thus implying that the bulk of merging 
activity for massive galaxies in clusters has to occur at much 
earlier epochs. 

This might appear not surprising as field studies have found 
evidence of a similar early assembly of massive galaxies with 
most of the stellar mass already assembled in systems more 
massive than t he loc a l character i stic m ass by redshift 1 (e.g., 
iFontana et alJ i20()A IConselicel §005)). Since galaxy evolu- 
tion in clusters is expected to be faster than in the field in hier- 
archical galaxy formation scenarios, an even milder evolution 
of the mass function of cluster galaxies is expected. The simi- 
larity between the shape of the MF we have found at z — 1 .2 in 
rich clusters and in the field sets upper limits on the difference 
of formation time-scales of stellar populations in high and low 
density environments. 

For example, if we use the star formation histories for early 
type galaxies in different environments, for different stellar 
masse s, as derived from the fossil record data (Th omas et alJ 
120051) . our results would imply an age difference of < 2Gyr. 



The situation is probably different for lower mass galaxies, 
however it remains difficult to probe the mass function signifi- 
cantly lower than 1O 1O M at z > 1. 

Among the caveats in our study, we should mention 
the so-called progenitor bias. Since we have considered the 
whole cluster galaxy population, our work is in principle 
not affected by the progenitor bias referred to when dealing 
with selected populations of galaxies (generally early-types 
(Ivan Dokkum & Franxl200ll) ). It remains true that galaxy pop- 
ulations in clusters at high re dshift mig ht not be directly com - 
parable to local ones (e.g., iKauffmann & Charlotl (Il998al) 1. 
and that the high-redshift clusters we are observing might 
not be the progen itors of the local X-ray luminous clusters 
dKauffmannll9 95). leading to an underestimated evolution. 

With the aid of multicolour photometry, including 
Spitzer/IRAC bands, we can now directly estimate the stellar 
masses of high redshift cluster galaxies, as well as approximate 
ages of their stellar populat ions, and push these studies out to 
Z = 1.4 ( lMullis etalJ(l2005l) . Stanford et al. (2005)), thus prob- 
ing an epoch which is thought to be crucial for the formation of 
massive clusters. This work will stimulate significant progress 
in discriminating between different formation scenarios. 
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